Towards an automatic ordering for a symmetric sparse direct solver
نویسندگان
چکیده
In recent years, nested dissection has grown in popularity as the method of choice for computing a pivot sequence for use with a sparse direct symmetric solver. This is particularly true for very large problems. For smaller problems, minimum degree based algorithms often produce orderings that lead to sparser matrix factors. Furthermore, minimum degree orderings are frequently significantly cheaper to compute than nested dissection. In this report, we look at whether we can predict which ordering will be better, using only the sparsity pattern of the matrix. Our aim is to choose efficiently a good ordering for a wide range of large problems from different application areas.
منابع مشابه
Effective Preconditioning through Ordering Interleaved with Incomplete Factorization
Consider the solution of a sparse linear system Ax = b when the matrix A is symmetric and positive definite. A typical iterative solver is obtained by using the method of Conjugate Gradients (CG) [15] preconditioned with an incomplete Cholesky (IC) factor L̂ [4]. The latter is an approximation to the (complete) Cholesky factor L, where A = LL . Consequently, the process of computing L̂ relies to ...
متن کاملPerformance of a Fully Parallel Sparse Solver
The performance of a fully parallel direct solver for large sparse symmetric positive deenite systems of linear equations is demonstrated. The solver is designed for distributed-memory, message-passing parallel computer systems. All phases of the computation, including symbolic processing as well as numeric factorization and triangular solution, are performed in parallel. A parallel Cartesian n...
متن کاملA Large-Grain Parallel Sparse System Solver
The eeciency of solving sparse linear systems on parallel processors and more complex multicluster architectures such as Cedar is greatly enhanced if relatively large grain computational tasks can be assigned to each cluster or processor. The ordering of a system into a bordered block upper triangular form facilitates a reasonable large-grain partitioning. A new algorithm which produces this fo...
متن کاملThe design and use of a sparse direct solver for skew symmetric matrices
We consider the direct solution of sparse skew symmetric matrices. We see that the pivoting strategies are similar, but simpler, to those used in the factorization of sparse symmetric indefinite matrices, and we briefly describe the algorithms used in a forthcoming direct code based on multifrontal techniques for the factorization of real skew symmetric matrices. We show how this factorization ...
متن کاملPSPASES: An Efficient and Scalable Parallel Sparse Direct Solver
Many problems in engineering and scienti c domains require solving large sparse systems of linear equations, as a computationally intensive step towards the nal solution. It has long been a challenge to develop e cient parallel formulations of sparse direct solvers due to several di erent complex steps involved in the process. In this paper, we describe PSPASES, one of the rst e cient, portable...
متن کامل